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Abstraction faite d'autres phénomènes importants, l'astronomie contemporaine doit surtout 
être caractérisée par le combat qu'elle livre contre les erreurs systématiques. 

Les erreurs systématiques sont ainsi appelées en opposition et avec les erreurs 
accidentelles et avec celles auxquelles on peut remédier en calculant des corrections, parce 
qu'on connaît la nature et le mode d'action des sources d'erreurs. Mais ces différences 
sont assez vagues. H y a des erreurs sur l'origine desquelles on peut avoir des doutes, 
encore longtemps après qu'on a trouvé la loi mathématique qui les représente comme une 
fonction de quelque circonstance se rapportant aux observations. Et si les limites entre 
les erreurs systématiques et les erreurs accidentelles sont incertaines, ce n'est pas seule- 
ment parce qu'il est douteux qu'on puisse, en thèse générale, attribuer des erreurs au 
hasard. En effet^ ce qui caractérise les erreurs accidentelles, c'est, d'une part, qu'en com- 
parant les résultats donnés par les répétitions d'une observation avec toutes les circon- 
stances qui l'accompagnent, on ne doit pas pouvoir constater de relation entre les 
erreurs et les circonstances en question, et, d'autre part, que, pour pouvoir être traitées 
comme accidentelles, il faut, lorsqu'on range les erreurs d'après leur grandeur et que 
le nombre des répétitions est suffisamment grand, qu'elles révèlent l'existence d'une loi 
qui mette la fréquence des erreurs tombant entre des limites arbitrairement choisies en 
relation avec ces limites elles-mêmes. 

Pour compléter la définition des erreurs accidentelles, on demandera peut-être que 
cette loi prenne la forme d'une fonction spéciale, à savoir l'exponentielle 



f(.x) == 






mV27r 



où ^0 ^^^ 1^ valeur la plus probable et m l'erreur moyenne. Mais même si l'on ne veut 
pas ranger parmi les erreurs systématiques celles dont la loi est autre que cette exponen- 
tielle, qui est la base d'un emploi direct de la méthode des moindres carrés, il y a des 
cas où la différence entre les deux espèces d'erreurs peut être très douteuse. 
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Dans plusieurs domaines, notamment dans celui des erreurs dites personnelles, qui 
semblent avoir une origine physiologique, et dans les séries d'observations servant à fixer 
les valeurs qui déterminent les particularités relatives à rinstallation et au mode d*action 
d'un instrument, il se présente des cas où Ton ne voit que trop clairement Timpossibilité 
d'exprimer la fréquence des différentes erreurs en fonction de leur grandeur, et où, lorsqu'on 
ordonne chronologiquement les observations de la même valeur, il semble qu'il se mani- 
feste une dépendance entre le temps de l'observation et les erreurs, en sorte qu'on pourrait 
être autorisé, en tout cas provisoirement, à appeler ces erreurs systématiques. Mais 
lorsqu'on a essayé de les corriger dans la supposition qu'elles étaient des fonctions mathé- 
matiques du temps de l'observation, l'espoir de débarrasser par ce moyen les observations 
des erreurs dont il s'agit a assez souvent été déçu. Lorsqu*on croyait avoir trouvé la 
formule de correction convenable et en faisait l'essai sur des séries d'observations conti- 
nuées pendant un temps suffisant, on aboutissait constamment à ce résultat qu'il était 
impossible de représenter les erreurs comme des fonctions du temps, ou du moins d'em- 
ployer la formule pour calculer à l'avance les erreurs personnelles ou les constantes des 
instruments. 

Peut-être qu'en pareil cas on obtiendrait souvent des résultats plus satisfaisants en 
déterminant plus exactement la relation de ces erreurs avec les erreurs accidentelles. Si 
les erreurs dont il s'agit ne sont pas simplement des erreurs accidentelles, il est au moins 
possible qu'elles résultent d'un concours de plusieurs sources d'erreurs accidentelles, car 
la proposition bien connue qu'une combinaison de différentes sources d'erreurs acciden- 
telles, agissant chacune selon sa loi, peut être représentée comme une source unique d'er- 
reurs accidentelles, avec une loi qui dépend de toutes les autres, présente une exception 
qui semble ne pas avoir été toi^ours suffisamment remarquée. Cette proposition cesse 
d'être vraie lorsque quelques-unes des sources d'erreurs exercent une action telle que, 
parmi les observations, qui sont d'ailleurs indépendantes les unes des autres, il y en a 
plus d'une qui est influencée par chaque erreur provenant de ces sources. Supposons, par 
exemple, qu'une série de valeurs o^ , o, , 03 . . . Oq, qui pourraient avoir été considérées 
comme Tobjet proprement dit des observations, aient, sous l'influence d'une source d'er- 
reurs accidentelles, reçu des accroissements e\ , e', . . . . e'. , d'après ime loi qui peut 
être supposée exponentielle, mais qu'on prenne pour objet de l'observation non les valeurs 
mêmes de mais les sommes successives de ces valeurs ordonnées suivant leurs indices, 
par conséquent 

Oj = Oi + o.^ 



et qu'en outre on soit forcé de procéder ainsi, parce que d'autres sources importantes 
d'erreurs accidentelles affectent ces sommes d'erreurs e'\^ «"^ , .... «"«, suivant une 
autre loi que nous pouvons également supposer exponentielle, les erreurs totales 



e*'^ "T" ^'i + ^'2 

(1) 

+ e\'\-e\+ ...«'„ 

ne suivront en général ni la loi exponentielle ni aucune autre loi, mais, lorsqu'elles seront 
ordonnées par rapport aux indices de e**^ présenteront des particularités qui ressemblent 
beaucoup à ce qu'on trouve dans les aperçus chronologiques mentionnés plus haut, con- 
cernant les erreurs personnelles et les résultats d'observation» faites sur les constantes des 
instruments, il y a aussi de très bonnes raisons pour expliquer les phénomènes en que- 
stion de cette manière ou d'une façon analogue par une complication d'erreurs acciden- 
telles, et, sous ce rapport, les erreurs personnelles se comportent à peu près comme celles 
qui proviennent des constantes des instruments, les organes des sens pouvant être consi- 
dérés comme des instruments. Les erreurs seront en partie des erreurs d'aperception, qui 
doivent varier indépendamment d'une observation à l'autre, de même que les erreurs e" de 
notre exemple, en partie des erreurs qui dépendent de l'état de l'instrument au moment 
de l'observation, et pour ces dernières il n'est pas en général permis de supposer qu'elles 
sont toiyours indépendantes de Tétat de l'instrument pendant les observations précédentes, 
en tant que cet état ne se laisse pas déterminer suivant les circonstances, comme la 
température, la pression barométrique, etc. 

Par contre, en ce qui concerne l'état de l'instrument, aussi longtemps qu'on se 
sert du même, il faut, à chaque observation, tenir compte de l'état où il se trouvait lors 
de l'observation précédente. On peut, il est vrai, pour beaucoup d'instruments, indiquer 
un état normal dont ils ne s'écartent jamais qu'en deçà de certaines limites, et, quand ils 
sont dérangés de leur position normale, il y a sans doute souvent plus de probabilité 
pour un retour vers celle-ci que pour un écart plus considérable. Cependant je suis très 
porté à croire que de telles variations accidentelles d'un instrument, qui le font osciller au- 
tour d'une position normale de manière qu'on doit toujours considérer celle-ci comme la 
plus probable (ces erreurs se combineraient alors avec les erreurs d'aperception «"), sont 
assez rares et, en tout cas, moins fréquentes que les variations où l'influence de l'état 
normal est insensible et où un «statu quo» est la règle, les variations accidentelles de 
l'instrument se produisant de manière que son état, à un moment donné, est aussi l'état 
h plus probable au moment suivant, d'où il suit que les erreurs se somment comme les e*. 
Quelques instruments présenteront peut-être encore une troisième espèce de variations 
qu'on pourrait appeler typiques à côté des précédentes, à savoir celles où l'état réel à 
chaque moment est égal à la moyenne entre l'état réel au moment précédent et l'état le 
plus probable au moment suivant. 

J'ai cherché à développer la théorie des erreurs quasi -systématiques qui, de la 
manière indiquée dans (1), proviennent de deux sources hétérogènes, et j'espère avoir con- 
tribué par là à réduire le domaine des erreurs systématiques. Assurément, il serait bien 
plus avantageux d'apprendre à calculer les corrections qui feraient complètement disparaître 
les erreurs systématiques que de les rapporter seulement à la catégorie des erreurs acci- 
dentelles, et d'indiquer des moyens pour l'emploi de la méthode des moindres carrés dans 



des cas où elle n'était auparavant pas applicable; mais ce n'est pas tout mal qui peut 
être extirpé avec la racine, et il est bien possible qu'on ait çà et là affaire à des erreurs 
dont les sources et les lois, pour être découvertes, exigent qu'on divise les erreurs ob- 
servées en deux classes, à savoir celles qui exercent leur influence sur les observations 
indépendamment les unes des autres et celles dont l'action s'accumule. 

Que cela soit possible, la preuve en est qu'il existe depuis longtemps une solution 
du problème pour un cas spécial. Si une série d'observations est divisée par de longs 
intervalles en un certain nombre de groupes, et si les observations dans chaque groupe se 
succèdent si rapidement qu'on puisse presque les regarder comme simultanées, le problème 
ne présente alors aucune difficulté; les répétitions faites dans les différents groupes servent 
à calculer, d'une part, les nombres moyens qui donnent lés valeurs les plus probables pour 
ces groupes, et, d^autre part, les erreurs moyennes qu'on doit s'attendre à rencontrer tant 
dans ces nombres moyens que dans chaque observation; quant aux valeurs les plus pro- 
bables tombant entre les intervalles des groupes, on les obtient par interpolation dans 
l'hypothèse d'une variation uniforme du groupé précédemment observé au suivant, et les 
différences qui se sont produites dans les longs intervalles peuvent toujours être utilisées 
pour juger de la stabilité de l'instrument. Le problème doit donc pouvoir aussi se résoudre 
lorsque les intervalles sont remplis par des observations éparses, el que les répétitions 
dans les divers groupes deviennent peu à peu de moins en moins nombreuses. La seule 
différence, c'est qu'il sera un peu plus difflcile d'obtenir les résultats désirés. 

Le principe pour la solution de ces questions n'est nullement artificiel: il faut 
seulement avoir soin, en mettant le problème en équation, de n'omettre aucune équation 
pouvant servir à compléter la description du problème, et de ne rien éliminer si ce n'est 
d'après les règles de la méthode des moindres carrés. 11 convient peut-être déjà de faire 
observer que le moyen analytique que j'emploie pour obtenir la solution, une fraction con- 
tinue, est bien un peu insolite dans les calculs de l'astronomie pratique, mais qu'il n'est 
en aucune façon incommode. 

Compensation d'une série d'observations d'une constante d'instru- 
ment. Nous supposons que x exprime quelque particularité d'un instrument (par ex- 
emple l'inclinaison d'un axe), laquelle serait constante si ce dernier était parfaitement 
stable, mais qu'en réalité l'instrument a varié d'un moment à l'autre, et que, par consé- 
quent, la «constante d'instrument», ^, au lieu de rester absolument constante, a subi des 
changements qui doivent être calculés par une série d'observations de œ. Nous suppo- 
sons en outre que chaque variation de œ^ d'un instant à l'autre, est purement accidentelle 
et suit une loi exponentielle déterminée par la condition que la valeur de sa est à chaque 
instant la valeur la plus probable pour l'instant suivant, tandis que la valeur moyenne des 
variations momentanées, /u, est connue. Si donc nous désignons par œ^^ œ^ , , . , Xn les 
valeurs que x doit être supposé avoir eues dans une série de moments <o ) *i . • • ■ ^« 
(en ordre chronologique), pendant lesquels on a fait des observations de œ ou pour les- 
quels, on désire de connaître x^ sera alors, d'après l'hypothèse, la valeur la plus pro- 
bable de A'rt+j — Xa et m^ = 2/z* (entre les limites a et a + M, le carré de l'écart moyen 

«.«+1 



pour cette différence, sera connu. Les poids absolus de ces équations auront pour va- 
leur V = -— ô- . On aura ainsi les n équations 
«,«+1 «* 

^1 — ^0 = 0> poids V 

0,1 

^2 — ^1 "^ ^> poids V 

*'^ ^ (2, 



•îtn — ^n-l= 0, poids V 

n— l,n 



Nous supposons de plus que les valeurs de a observées dans les temps ^q? ^i • • * ^n 
sont affectées d'erreurs accidentelles, suivant des lois exponentielles, et que nous con- 
naissons les valeurs les plus probables (avant la compensation) ^^ , z^ , , . . z» prove- 
nant des observations, ainsi que les erreurs moyennes m|, m^ . . . . m» ou les poids t;^, 

î?2 . . • . v«, en général Va = — r» On obtient ainsi les n-f 1 éc^uations 

^0 — ^0 = ^> poids Vq 
z — ^. = 0, poids v, 

. (3) 

Zn — Xn=^ 0, poids Vn 

qui, conjointement avec (2), renferment toutes les conditions Indépendantes les unes des 
autres qui doivent servir à déterminer les n-\~\ inconnues Xq^ œ^ •••^n, que nous nous 
proposons de trouver avec leurs erreurs moyennes après la compensation. 

Remarque. SI, parmi les temps t, il y eu a quelques-uns où Ton n*a pas fait d'observation, et 
pour lesquels on demande seulement de connaître les valeurs les plus probables de xa. il faudra poser les 
poids correspondants va = 0, ce qui rend les valeurs de za arbitraires (mais non infinies). 

Si maintenant on traite les équations (2) et (3) par la méthode des moindres carrés, 
il en résulte les n-[-\ équations : 

0,1 0,1 

"^l^l "^^ Vd7oH~ (^"|~ ^1 + *^)^1 — ^^2 

0.1 0,1 1,2 1,2 

— VXa-X^(v + Va'\-V) Xa V Xa±i > (^) 

a — 1,« a — l,a «,«+1 a,a4"l 



VaZ. 



a^a 



Vn^iZn-i= —V A'„_2+ (« + Vn-1 + ^) ^«-1— W ^« 

n— 2,11 — 1 n — 2, n — 1 n — l.n n — l,n 

n— 1,» n — l,n 

Dans ces équations il faut d'abord éliminer les x» correspondant à Va = 0, par conséquent 
aux observations fictives, qui n'y figurent que pour servir à la détermination de la valeur 
de X eii dehors des temps des observations réelles. Mais de t'a =-^ il suit que 
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(ara — «a-i):— = |ir«+i — «,) : — = («.+i — ««) : (— ■ + — ) <S) 

V V \^ ^ / 

a— l,fl tftfl+l a— l.a a,«+i 

de sorte que Xa pourra être déterminé aussitôt qu'on aura trouvé ^Ta-i et œ^^ par une 
simple interpolation entre ces valeurs, dans la supposition que les variations de x sont 
inversement proportionnelles aux poids entre les limites de Tintervalle considéré, ou directe- 
ment proportionnelles aux carrés des écarts moyens ou à Tintervalle de temps lui-môme, 
en tant qu'il y a eu à chaque moment une chance égale pour un déplacement de Tinstru- 
ment. Si la première ou la dernière observation, Zq ou z^^ est fictive, et par conséquent 
si Ton demande les valeurs de la constante de Tinstrument avant et après la période où 
il a été fait des observations, il suit alors de t;o ==- que ^o ="^i> ^« = ^î ^» = ^n-i^ et 
les valeurs de œ qu'on obtient pour le premier et le dernier moment des observations réelles 
sont donc les valeurs les plus probables respectivement pour les temps qui précèdent et 
pour ceux qui suivent. 

Les équations (4) montrent en outre, ce qui aussi va de soi, que Télimination 
de œay lorsque fa = 0, n'apporte dans ces équations d'autre changement que de réduire 
le système des équations à ce qifon aurait obtenu immédiatement en omettant l'observation 
fictive; on n'avait qu'à poser 

m^ =^ m^ -{- m* 

a— l,a+l 41— 1, a 9ta-\^ 



OU V 



V V 
a— l,a a, a-)-l 



o— 1,0 a,a+l 

et à passer directement de l'équation 

a— 2,a— 1 «— 2,«— 1 a—Un-^i u — li«i-l 

^ Vu-^iZa-f-i = —V Xa^i-\-(V+Va^i+V) .r^+i — t? Xa^2 

a~-l,a+l a-l.a+l fl+i,fl+2 «+-1,0+2 

en laissant de côté l'équation de VaZ^. Nous pouvons donc, lorsqu'il s'agit de calculer de 
véritables exemples, nous dispenser d'introduire des observations fictives, et nous savons 
maintenant comment, avec les valeurs de x correspondant aux temps des observations réelles, 
on peut trouver x pour un autre temps quelconque. 

Les équations (4) doivent naturellement être résolues de la même manière 
que les équations normales de la méthode des moindres carrés, et ce travail est faci- 
lité ici par la forme relativement simple des équations; mais, pour en tirer parti et pré- 
senter tout de suite la solution sous la forme la plus commode pour le calcul, il faut 
avoir recours à quelques quantités auxiliaires. Nous devons d'abord, au lieu des poids 
Va et Va,a+i, întroduirc deux autres séries de nombres, déterminés par les formules récur- 
rentes : 

Wa+l = Vai-i + U (6) 

a,«-hl 

et ± = ±_|_± (7) 

U V u„ 

a, â-jri «i«-hi 
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En supposant Uq = t?o > ^^ calcule facilement Ua et u jusqu'à u et Un à l'aide d'une table 

a,a-f-i n— l,n 

de nombres réciproques. 

On déduit ensuite des z une série de valeurs yo ) ^i • * • y») Qui, dans la suite du 
calcul, doivent remplacer les z^ et dont chacune représentera en même temps la valeur 
correspondante de /u, dans Thypothèse que toutes les observations avec des indices plus 
élevés ont été laissées de côté, d'où, en tout cas, y» = ^n* Ces valeurs sont données 
(pour des indices croissants) par la formule récurrente 

ou plus commodément par 

Wa+l (y«+l — !/a) = Va+i (Za H — I/a) ] 
OU J (8) 

«r,a+l 

conjointement avec yo = -^rQ. 

En opérant ces substitutions dans (4), on trouvera que les valeurs de a sont déter- 
minées par la formule récurrente (pour des indices décroissants) 

W«ya = (î/a + V) ^« — V aa^i (9) 

OU plus commodément par 

V (Xa — ^a+i) = M (ya — «a+i) (9 ») 

OU «•«+! fl»a+l 

Wa(^a— yj = w (.r«+i — yj l9b) 

conjointement avec a^. = y«. 

Que non-seulement cela soit exact, mais que (9) soit justement aussi la forme type 
du résultat de l'élimination suivant la méthode des moindres carrés, c'est ce qu'on peut 
démontrer par la méthode d'induction générale. Car, pour a = 0, (9) est identique 
avec la première équation (4), et si l'on suppose que (9) est valable pour tous les indices 
jusqu'à a, et, par conséquent, que l'élimination de jtq, ^i . . .^«-i a conduit à un système 
de « — a+t équations, à savoir (9) et les n -- a dernières (4), l'élimination de ^a entre 
les deux premières de ces équations 

Va+l^a+i= — V ^« + (u + Va-\-i + V) IVa fl — ^ ^«+2 , 

les seules qui renferment encore ara, se fera en multipliant la première par 



V u 

a, ii+l g, o-hl 

a. «+1 



(voir (7)) 



et en ajoutant le produit à la seconde, ce qui donne d'après (7) 

W yfl+ Va+iZa^i = (m + î^a-t-l't" ^) ^^fi ^ ^ ^«"1-2 
a,o+l a.a+1 a-f-1,a-h2 €i+l,«+2 



10 

et d'après (8) et (6) 

«+l,a+2 «-M.IC+2 

dont raccord avec (9) prouve rexaclitude de la proposition, à condition seulement que 
a<n. En outre, on voit facilement que la dernière élimination (de «««i) conduit à Unt/n 

En développant l'expression explicite de m„ au moyen de (6) et de (7), on trouve la 
fraction continue 

,1 



n— i,» 






tJij Vl+ ï 



0,1 



et, si Ton rejette un nombre quelconque des premiers termes, les fractions continues for- 
mées par les termes restants seront alternativement u„ et — • 

Vu l'importance que cette fraction continue a pour notre problème, on peut dire 
avec certitude que des solutions plus simples ne deviendront possibles que lorsque la frac- 
tion continue prendra une forme plus simple, et par conséquent qu'on obtiendra des solu- 
tions approximatives d'un calcul plus facile seulement dans les cas qui s'en rapprochent. 
Comme remplissant cette condition, on pourrait d'abord penser au cas où la fraction con- 
tinue devient périodique, mais il n'en résulte aucune simpliflcation sensible dans le calcul 
pratique. Il en est tout autrement dans le cas où la fraction continue est interrompue 

parce qu'un dénominateur partiel devient infini, c'est-à-dire lorsque »« = oo ou r = 0, et 

«.«+1 

dans celui où la fraction continue est réduite à un nombre plus petit de termes parce 
qu'un dénominateur partiel devient nul, c'est-à-dire lorsque r. = ou r = oo . 

Lorsqu'une observation peut être considérée comme exemple d'erreur, ©« = 00, ou 

que la série des observations est interrompue par un long intervalle ou un intervalle 

pendant lequel on peut craindre que l'instrument n'ait subi un dérangement, i; = 0, la 

«.«+1 

fraction continue est interrompue ou se divise en deux fractions continues indépendantes 
Tune de l'autre. Il en résulte que le calcul des a se fractionne aussi en deux parties; 
dans le premier cas, de manière que u ^= v et Za ^^ 2/a == ^a, et, par suite, que les deux 

«»a+l ff,«+l 

séries sont reliées entre elles par la valeur commune Xa = Za\ dans le second cas, de 
manière que Wa+i = r^+i et x^^^^ya^ mais ya^i = *«+!) en sorte que les deux séries 
n'ont entre elles aucune liaison, et sont séparées au moins par une valeur intermédiaire 
de /c, qu'on ne peut déterminer. Des firactionnements de ce genre pourront dans beaucoup 



Il 



d^autres cas avoir un grand avantage pour le calcul; mais celui qu^on en retire ici n'est 
pas grand, car on pourra presque toujours fractionner arbitrairement le calcul sans autre 
inconvénient que celui-ci, qu*un double calcul peu compliqué des u, des y et des x autour 
du point de partage aboutira rapidement à relier les parties qui avaient été traitées séparé- 
ment. Ces considérations prouvent seulement que, lorsqu'on a le choix, il faut fractionner 
de manière à faire entrer dans les deux parties une valeur bien déterminée de z^ ou inter- 
rompre le calcul en un point où il y a un grand intervalle entre deux observations consé- 
cutives, dont Tune sert à terminer la première partie, et Tautre à commencer la seconde. 

Relativement aux cas où le nombre des termes de la fraction continue est suscep- 
tible de réduction, nous pouvons nous limiter à celui de v = oo , comme nous avons déjà, 

tf.a+l 

à un autre point de vue, considéré le cas Va = et vu que, là où il s'agit de trouver osa 
pour une époque pour laquelle il n'existe pas d'observations (ou seulement des observa- 
tions méritant peu de confiance et n'ayant qu'un poids insignifiant), on doit calculer les 
valeurs les plus probables au moyen d'une interpolation simple entre les deux époques 
d'observation les plus voisines qui ont précédé et suivi l'époque en question, dans l'hypo- 
thèse d'une variation uniforme dans l'instrument. Par contre si v = oo , ce qui suppose 

que les temps t^ et <a+i des observations se succèdent si rapidement que la position de 
l'instrument ne peut avoir varié dans l'intervalle, il est évident qu'on est autorisé, en ayant 
égard aux poids, à former des moyennes de ces observations consécutives et, par conséquent, 
à considérer 



+ *^«-H 



avec le poids t?a + Vi,+i 



comme observé au lieu de Za et Za+i ) ce qui est aussi confirmé par notre théorie, laquelle 
donne pour t? = oo 

tf— l,a 
a— 1,0 

et *ptf+i = ^. 



a • 



Si, dans ce qui précède, j'ai eu soin de ne pas m'écarter du procédé d'élimination 
propre à la méthode des moindres carrés, c'a été pour pouvoir tirer avantage, dans les 
déterminations des erreurs moyennes, de l'importante proposition que je connais seulement 
par les communications verbales de M. le professeur Oppermann (congrès des naturalistes 
à Copenhague, 1873), à savoir, d'une part, qu'il existe des systèmes de fonctions linéaires 
d'observations qui, en général, peuvent être considérées comme des observations indépen- 
dantes les unes des autres, et, comme telles, prendre la place des observations réelles, 
notamment dans tous les calculs concernant les erreurs moyennes des fonctions des élé- 
ments; et, d'autre part, que les premiers membres des équations résultant de l'élimination 
par la méthode des moindres carrés, ici les équations (9), représentent précisément un 
pareil système de fonctions qui peuvent être prises pour des observations indépendantes les 
unes des autres. 

2* 
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Par conséquent, dans les équations 

0,1 0,1 

1.2 1,2 

M„_iyn-i= (w„_i+v) a?„_i — t? a?„ 

»— l,n n— l,it 

«n yn = Wi.^n 

les produits uq^o j • • • ^a^a > • • • ^nyn sont à considérer comme des valeurs observées indé- 
pendamment les unes des autres, de manière que le carré de Terreur moyenne de Uaya est 

^(^uy*) = Ua-\-v 

a,a+l 

par suite les valeurs de y elles-mêmes peuvent aussi être regardées comme des observa- 
tions indépendantes les unes des autres, de sorte que Ton a pour les carrés de leurs 
erreurs moyennes 

^(y«)=— ^ = — ' — (H) 

avec la seule exception que 

M^iyn) = APM =— (Ha) 

^» 

Si maintenant on repréîsente explicitement les a comme des fonctions linéaires des y, il 

sera très facile de calculer Terreur moyenne d'une valeur quelconque provenant de la 

compensation, et donnée d'avance comme fonction linéaire des a. Mais comme 






les coefficients des équations 

^0 = cyo+ cyi + cy2+ . . , + c yn^i + cyn 

0,0 0,1 0,2 0,fi— 1 0,n 

A'i = cyi-\-cy2-\-...-^c yn-^i+cyn 

1,1 1,2 l,fi— 1 l,n 

^2 = <?y2+...+ C yn-l+Cy« 

2,2 2,11^1 2,1» 



(12) 



^11-1= cyn-^i + ^y» 

n— 1,11 — 1 n— 1,« 



3:n = y. 

seront déterminés par 

u 

U U 

<i,«+l «+l,a+2 6—1,6 



C = 

et par c = pour i<a. 

«,6 



±i(l^^tl) = c ^--+îpour6>a 

j_l \ Uô / a+1,6 W« 



(13) 
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Pour leb carrés des erreurs moyennes on aura donc 

et en général 

— « ffft"* — ** V^-o «•*/ 

de la première de ces formules on déduit la formule récurrente (pour des indices décrois- 
sants) 

^^"'^ = -^ + [!^) m^^,) (14) 

Comme cas spécial de la formule générale, nous remarquons surtout celui qui donne le 
carré de Terreur moyenne d'une fonction linéaire rxa + «^a+i de deux a consécutifs 

Ua — u / ^' \^ 



en effet, comme fonction linéaire seulement des y avec des indices plus élevés, a^^i peut 
être considéré comme une observation indépendante à côté de ya* En posant r »» 1 et 
8 = — 1 , il vient 



Ua — u /Ua — u \2 



(15) 



Les équations (14) et (15) sont très commodes pour calculer les erreurs moyennes qui, 
après la compensation, peuvent rester dans les équations (2) et (3), lesquelles résultent 
directement des conditions fondamentales du problème. Avec leur aide on peut aussi 
transformer Téquation générale de ]kP{raa-{-8a!a+i) en la suivante 

AP(rXa + «^«-h) = (»• + «) {r^M + 8M^(Xa^i)) — r8AP(Xa — ^«+i). 

Si maintenant on se rappelle que, suivant (5), les a pour lesquels il n*existe pas 
d'observation réelle sont déterminés par les x qui correspondent aux observations précé- 
dente et suivante les plus voisines, comme étant une fonction linéaire de ces derniers — 
si Ton veut rxa + «^a+i , '^ + * ^^Qt égal à 1 — on sera très porté à employer cette 
valeur de J/^ (r^^a + ^^«-h) dans le calcul de Tincertitude moyenne des x qui ne sont 
pas déterminés par des observations correspondantes. Mais ce serait une faute grave dont 
il faut bien se garder dans des recherches de ce genre. Car même si Xf, «= rxa -[- 8Xa^ , 
il n'est pas permis pour cela de poser M^(x^) = iP(rxa-\'8Xa-\-i) dans les cas où l'équa- 
tion x^ = rxu -f 8Xa-^i est elle-même le résultat de la compensation et par suite n'a pas 
une valeur absolue, mais est seulement valable pour une approximation ayant pour mesure 
une erreur moyenne dont le carré sera M^(x^) — Jf * (r^'a + «^a+i). M^(Xb\ doit donc être 
déterminé directement par le même moyen qui a donné «», à savoir l'introduction d'une 
observation fictive. Les conséquences de cette introduction s'expliquent peut-être le mieux 
par le fait que les schémas primitifs des v, des u, des y et des x 
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v^ 



V 

tf,tf+l 



Vai-i 



U. 



U 



Ua^i 



y 



a 



se changent en 



V, 



V 



V6 



V 



u 



a 



U 



Ut, 



u 



Va+1 



conformément aux relations suivantes 



Wa+l 



ya+i 



y« 



^« 



a?. 



a-|-l 



tf* 



y«+i 



iTx 



Ot?» 



•îPa+l 



V* = 

1 + 1 = ± 



Ut = U 



UaV 
a, 6 



V U 
*,a+l a,a-f-l 






«1,6 «a+V V — W 

a, 6 6,a+i o,«-|-4 

- U 



y*=y« 



œb = 



a, 6 6,0^-1 



OÙ par conséquent 



P 



V 



a,6 M+1 «.«+1 



et 



l-p = 



Î7 

r-f- î7 

a, 6 6,<ff-l 



U (w« — Wé) 
a,al-l 



Relativement aux y considérés comme des quantités observées indépendamment 
les unes des autres, tout le changement se réduit donc à ce que Vifa primitif, qui avait le 

poids Wfl — M, est remplacé par deux autres, à savoir i/a avec le poids Ua — u^ et y* =y« 

«,«+1 

avec le poids u^ — u, de sorte que les valeurs elles-mêmes sont identiques et la somme des 

poids est égale au poids primitif. Mais, quelque petit que paraisse être ce changement, 
il a cependant une signification réelle, mais seulement par rapport à la détermination de 
AP(Xb), En effet, d'après (14), on a directement 
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«4 U / ** \ * 



OD a en outre par un double emploi de (14) 






qui, pour ti» = u, donne 



Ua — U / " V 



c'est-à-dire Téquation (14). On a donc 

pAfiM + a -p)MHa!,^i) = — ^^ + '-^ '-^ m^a^i) 

et 

(Ua — u^) (w* — u) 

JiP(x,) = pilf»(a!.) + (1 -p)AP{a^.^i) + ^-î^^ (1 - u 3f»(^.+ir) 

(u,-« )»/ 1 \ 

par conséquent, suivant (15) 

M^(xt) = pilf«(«.) + (1 —p)Mi{x,+i) +/>(! -/>) ( — î!^ — A/«(^. - ar„^.t) ) 

ce qui donne pour résultat 

Si, dans ces formules, on fait t? = 0, cas dans lequel on aura à déterminer les a 

et leurs erreurs moyennes en dehors de la période que comprennent les observations, où 
dans un intervalle où Tinstrument peut avoir subi une variation telle que la série des ob- 
servations doive être considérée comme complètement interrompue, il en résultera suivant 
(7) qu'on devra aussi avoir u = 0, Si Ton n'a pas en même temps v == 0, on aura 

p = \^ xt = Xa et M^(xt) == — = 1 , ce qui est d'ailleurs naturel. 



H 



CumcnuCaBL j» ok voÊBgtiBKSiaM. juiaiiiiiia. 
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TÇT 
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Co — 1^0 = Co — oq et on place ensuite ^^ au-dessous de Oit à une distance égale à celle dont rjo précède Cn 

1 1 Ui 

ce qui donne «i — ç?i = d — j^q c^» par conséquent, diaprés (6), «i — Ci = 1 = • SI Ton mène 

0,1 ^ 0,1 

alors la droite i^of'i ^^ parallèlement à celle-ci par Vi, la ligne t^iU^, ^7, est déterminé de manière que 

Cl — î7i = 1 :«i et Ca — yi = I :w. 

1,2 

En continuant ainsi, aussitôt que r^a est trouvé, on construit ^a+i au-dessous de e/a+i de façon que 

0a4-i — ç>a+i = C«+i — î?«= 1 :t>+l:tta = l:tt (suivant 7), 

et par suite t/a+i — Ca+i = — — (suivant 6); 

puis on construit l'abscisse ija-\-i en menant ^^-i-if/a+i =6:7«^a+if et ainsi de suite jusqu'à ce que r^n soit 
trouvé (sur la figure la ligne ijiUi coïncide par hasard avec une partie de 7i?>a)* 

Cj Après avoir élevé des ordonnées verticalement au-dessus de tous les points t^», on pourra y 
construire les points ya , qui , par la longueur ya — rjtiy. représentent les quantités auxiliaires ya . On cons- 
truit d'abord ^0 de façon que ^q — 7o = ^^q — Coi et Joint ensuite 2^0 ^ ^t P&r une droite qui coupe précisé- 
ment l'ordonnée de t^^ en y^ . Sur la droite j/jS, est situé y, et ainsi de suite; en général y^-^-i se trouve 
à l'Intersection de l'ordonnée de rja-^-i avec la droite qui joint y a et za+i (voir les équations 8). 

Lorsque les points ya ont été déterminés jusqu'à yn , on peut, à l'aide d'une construction rétrograde 
successive, représenter les résultats de la compensation par une série de points xa, qui, étant situés verti- 
calement au-dessus ou au-dessous des points za, c*est-à-dlre comme ceux-ci sur l'ordonnée de Ca» peuvent être 
considérés comme les positions compensées des points 2*. Parallèlement à Taxe des abscisses on mène par 
yn une droite qui coupe l'ordonnée de Cn au point Xn] puis on Joint les points yn^i et Xn par une droite 
dont l'Intersection avec l'ordonnée Cit— is^n—i donne le polnt.arn— 1> et ainsi de suite, la ligne qui joint y» et 
Xé-^-i rencontrant généralement l'ordonnée C«2a au point xa» conformément aux équations (9 a) et (9 b). Le 
dernier x, c'est-à-dire Xq, une fois trouvé, on peut encore continuer en menant par a^o ^^^ parallèle à l'axe 
des abscisses. 

Ce n'est pas seulement aux points Xq, Xi ... xn-^i, mais aussi, suivant (5), pour d'autres abscisses 
quelconques et par conséquent pour un temps quelconque, que la ligne brisée a^o» ^n ^3 * • ^ ^n, addition- 
née des lignes parallèles à l'axe des abscisses pour les parties de — 00 jusqu'à Xq et de xn à l'Infini, repré- 
sentera les variations de la constante de l'Instrument de Ja manière qui, d'après la compensation, doit être 
regardée comme la plus probable. SI les observations ont à l'origine été rendues positives par l'addition 
d'une constante arbitraire, tl faudra naturellement la retrancher des ordonnées Xa — Ca lorsqu'on voudra me- 
surer la quantité xa. 

Détermination des unités pour les deux espèces de poids. 

Après n'importe quelle compensation, surgit toujours la question de savoir si les 
erreurs qu'elle laisse subsister dans les observations satisfont aux suppositions qu'on a 
faites sur leur nature, et sur lesquelles la compensation repose; si Ton ne peut ou ne 
veut pas la résoudre, le résultat ne saurait être considéré que comme hypothétique.- H 
faut entre autres chercher à éclaircir si les suppositions, relatives à des erreurs moyennes 
ou à des poids, sur lesquelles est basée une compensation d'après la méthode des moindres 
carrés, s'accordent avec ce qu'on peut déduire à posteriori du tableau des erreurs restantes. 
La question se présente souvent sous cette forme particulière, qu'on ne connaît pas à priori 
les poids des différentes équations de condition, mais bien quelques systèmes de rapports 
existant entre eux, desquels seraient sufûsants pour la détermination des poids, si l'on 
connaissait en outre un petit nombre d'unités de poids, une pour chaque système; il faut 

s 



18 

donc ranger ces unités de poids parmi les inconnues du problème, et, par suite de la 
nature compliquée de ce dernier, on en est réduit à la méthode indirecte de commencer 
avec des unités de poids purement hypothétiques, et, après s*en être servi pour calculer 
la compensation, de déduire des erreurs restantes des valeurs nouvelles et en général plus 
approchées des unités de poids, valeurs avec lesquelles on recommence le calcul jusqu*à 
ce que les unités de poids à posteriori après la dernière compensation présentent un accord 
suffisant avec les unités de poids à priori avant cette compensation. On n*échappe entière- 
ment à ces répétitions que dans un seul cas, à savoir lorsqu*une unité de poids commune 
à toutes les observations a été inconnue; dans ce cas, la valeur absolue de Tunité de 
poids est sans influence sur le résultat de la compensation, et on la détermine par la 

formule 

n — m 

où / est Terreur restante dans une des équations, n le nombre de ces dernières et m 
celui des éléments, ou peut-être par 

n — m — 1 

si Ton doit, sous tous les rapports, ranger Tunité de poids inconnue parmi les inconnues 
du problème, ce qui cependant est un peu problématique. Les exemples de cas avec plu- 
sieurs unités de poids inconnues sont très fréquents en astronomie; il suffit, à cet égard, 
de nommer les poids des ascensions droites et des déclinaisons dans la détermination 
d'une orbite. Les exemples d*un traitement exact de cas avec plusieurs unités de poids 
inconnues sont assez rares. 

Dans les problèmes spéciaux que nous traitons ici, on sera ordinairement forcé de 
recommencer au moins une fois la compensation avec des poids plus exacts. Car, tandis 
qu*on connaîtra assez souvent les rapports entre 

et entre 



Vqj ro, 


... Vn 


t;, t>, 

0,1 1,2 


. . . V 

Il — 1,« 



on ne connaîtra pas toujours assez exactement Tunité des poids v^,, et le plus souvent on 
sera placé dans des conditions encore plus défavorables vis-à-vis de Tunité des poids v . 

«»a+l 

Par suite notamment d'une circonstance que nous allons mentionner, on ne pourra éviter 
de recommencer la compensation; cependant il semble y avoir lieu de supposer que la 
convergence dans la détermination des unités de poids sera en général si grande que le 
travail ne deviendra pas facilement impraticable. 

Lorsqu'une compensation a été calculée sur la base d'une unité hypothétique E^ 
pour les poids t?o, vi . . . v«, et d'une autre E„ pour le groupe t?, v, . . . r, il est évident 

0,1 1,2 fi— 1,11 

qu'il faut corriger E^ par les difi'érences Za — ^a et E„ par les différences Xa — ^«+i. 
Nous avons n+ 1 différences de la première espèce, n de la seconde, et n-f 1 éléments 
^0, ^1 ...^Pi sont déterminés par toutes les 2n+l équations. A ces éléments il faut 
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peut-être encore en ajouter 2 (pour E^ et E^)\ mais, à cause des grandes valeurs que n 
aura toujours dans des recherches de ce genre, cette question devient tout à fait secon- 
daire, et je ferai dans ce qui suit abstraction de cette possibilité en m*en tenant à n+1 
pour le nombre des éléments. On peut donc pour les unités de poids cherchées poser 

j^ ^ n+ 1 — f jg7 ^ n — 7i 

«,«+1 

OÙ cependant on ne sait encore autre chose de f et de g;, si ce n*est que leur somme 
représente le nombre des éléments 

^+r) = n+ 1; 

il s'agit de déterminer quelle est la partie de ce nombre qui appartient à chacun des deux 
groupes. Dans d'autres problèmes où les équations sont très nombreuses relativement 
aux éléments, cela n'a pas grand* chose à dire, et on pourra souvent se contenter de 
répartir également le nombre des éléments; mais ici on n'a jamais deux fois plus d'équa- 
tions que d'éléments, et un exemple montrera facilement que la répartition doit quelquefois 
être très inégale. Supposons en effet qu'une constante d'instrument n'a été observée qu'à 
deux époques essentiellement différentes, mais chaque fois avec de fréquentes répétitions 
se succédant sans interruption et sans qu'on ait eu à craindre un changement réel dans 
l'instrument; il sera alors facile de voir que tj se rapprochera de la limite n — 1, tandis 
que f se rapproche de 2. En général, il sera donc nécessaire, pour déterminer f et 3; 
dans des cas plus compliqués, d'examiner en détail, pour chacune des équations de condition, 
dans quelle mesure elle contribue au nombre f ou 37, qui doit être retranché de celui des 

équations, ou combien il faut ajouter respectivement à Va(Za — ^a)* ou à v (^« — a?fl+i)* 

«.•+1 

pour avoir le droit de supposer que la somme obtenue doit en moyenne être égale 4 1 
pour chaque équation. Sous la dernière forme de la question, la réponse est facile; car 
après les compensations , Xa et Xa — Xa^i sont aussi affectés d'une incertitude ayant 
respectivement pour mesure les erreurs moyennes M(Xa) et M(Xa — a?«+i), voir (14) et (15). 
Comme cette incertitude, par suite de la circonstance que la somme des carrés des écarts 
multipliée par le poids doit par la compensation être réduite à son minimum, est juste- 
ment cause que cette somme de carrés est moindre que si les observations avaient été 
comparées avec les valeurs vraies, il est du moins très vraisemblable qu'on opère exacte- 
ment en ajoutant au carré de chaque écart le carré de Terreur moyenne qui -reste après la 
compensation. C'est en procédant d'une manière analogue qu'on démontre sommairement 
la règle générale, que le nombre des observations doit être diminué de celui des éléments 
dans la formule qui sert à déterminer a posteriori les poids. Dans le cas qui nous occupe, 
on verra sans grande difScuIté à l'aide des formules (12) et (13) que 

IVaM^(Xa) + 2V iP(Xa— ^«+l) = n + 1. 

Pour les deux groupes d'équations on trouve de cette manière 

f = 2v^M^(x^) 



Tj = 2v iP(Xa—Xa-Yi) 

«fO+1 



3* 



_^0 

et en remarquant que ces valeurs de f et de ^ sont indépendantes de chaque multiplicateur 
commun de E^ et de E„^ et ne varient ordinairement que peu par une variation de E^:E^^ 
nous aurons pour la détermination des unités de poids 

n — 2v idP(œa — A'fl+i) ) (17) 
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Le signe indiquant que nous avons calculé avec des poids justes et absolus la dernière 
répétition de la compensation sera donc E^ = \ et J&» = 1 . 

Il faut remarquer ici que, tandis que Va iP(Xa) ne semble pas pouvoir être simpliûé 
avec quelque avantage par (6) et (14), on a suivant (7) et (15) 

équation qui nous fait voir que, lorsque v est relativement grand et par suite »« — ^ 

petit, V M^(Xa — ^a+i) sera en général un peu plus petit que 1, en sorte que les différences 
«,«+1 

Xa — ^a+1) ^^^^ l6s cas OÙ Ics tcmps ta ct ta^t dcs obscrvatious diffèrent peu les uns des 
autres, ne pourront pour ainsi dire contribuer en rien à la détermination de E^. Il suit 
de là que les valeurs compensées de^a et de Xa^i^ dans des cas pareils, ne présenteront 
guère de différence sensible, et nous en tirons la règle générale qu'avant de procéder à 
un.e compensation de ce genre, on peut et on doit réduire tous les groupes d'observations 
de cette espèce à de simples moyennes séparées seulement par de petits intervalles; avec 
un peu de précaution on évitera facilement de pousser trop loin cette réduction, qui évi- 
demment épargne beaucoup de travail. Naturellement, dans la détermination de Tunité de 
poids J?r, il ne faut pas négliger d'utiliser les observations qu'on réunit dans les moyennes. 
En opérant ainsi, on réussira souvent, avant la compensation proprement dite, à obtenir 
une si bonne détermination de E^ que la compensation finale ne lui fera subir qu'une 
modification insignifiante. Il est plus difficile de s'assurer à l'avance d'une détermination 
suffisamment approchée de ^». A cet égard, je conseillerai de faire un calcul provisoire 
en prenant pour principe de pousser d'abord jusqu'à Textréme limite permise la réduction 
d^s observations consécutives en moyennes, et de considérer ensuite ces moyennes comme 
des valeurs compensées de x^ en faisant complètement abstraction des observations éparses 
qui ne peuvent servir au calcul de moyennes ; on se procure ainsi une série de valeurs 

Xa^ Xb^ Xc^ XDy Xe* 

En même temps il faut réunir les poids v fixés provisoirement avec une unité arbitraire, 

c.a-f-l 

de manière que 

1 ^1 

— = y — y etc. 

A,B «.«+1 



;» •--*-^i 
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Si l'on calcule alors l'unité de poids E^ d'après la formule 

iV+ 1 étant le nombre des moyennes xa employées etc., cette détermination sera en 
générai un peu trop grande, mais pourra cependant être utilisée comme une première 
approximation. 

Mais, même si Ton a été conflrmé par une compensation dans Tidée que les poids 
employés sont justes, dans Thypothèse que les lois des erreurs sont des exponentielles, 
il reste encore pour la critique plus de place qu'on ne le croit généralement. Les lois 
asymétriques ou les lois avec plusieurs maxima se présentent, il est vrai, rarement; mais, 
le cas échéant, elles exigent une autre manière de procéder. Quant aux erreurs systéma- 
tiques, il n'est pas facile de s'en délivrer complètement. 

, Compensation d'après des fonctions linéaires d'éléments inconnus. 

Lorsque, en suivant les règles exposées plus haut, on aura terminé une compensa- 
tion pour une série d'observations avec complication de deux sources d'erreurs acciden- 
telles de nature différente, on ne sera naturellement pas dispensé de soumettre les deux 
tableaux des erreurs ainsi trouvées à un nouveau contrôle. Il est en effet possible que 
les résultats présentent des manques indiquant que la complication est d'une autre espèce 
que celle relativement simple dont il est question ici. 11 se peut également qu'après avoir 
séparé les deux espèces d'erreurs, on soit en état de découvrir dans les tableaux des 
traces d'erreurs systématiques proprement dites, erreurs qu'on ne pouvait constater tant 
qu'elles étaient cachées par la complication entre les erreurs d'observation proprement 
dites et les déplacements de l'instrument, mais qui apparaissent ensuite si distinctement, 
qu'on peut avoir l'espoir de les éliminer en calculant des corrections basées sur quelque 
circonstance secondaire se rapportant aux observations. D'ailleurs nous avons jusqu'ici 
seulement considéré le cas où des observations chargées d'erreurs compliquées concer- 
naient un phénomène qui en soi pouvait être supposé constant. Il nous reste donc encore 
à exposer en peu de mots comment on doit s'y prendre, lorsque le phénomène observé 
varie suivant des lois qui peuvent être exprimées par des fonctions linéaires avec un certain 
nombre de coefQcients inconnus, en même temps que les erreurs sont de l'espèce com- 
pliquée dont il est question ici. Il suffira de traiter le cas où la loi du phénomène ob- 
servé est Za = pfa-^qffa, équatiou dans laquelle C«) qui est la valeur vraie de l'observa- 
tion faite dans le temps ^«j dépend de deux éléments inconnus p. et j, tandis que /« et 
ffa sont des nombres connus qui dépendent de ta. Les observations sont supposées affec- 
tées de deux espèces d'erreurs, les unes, Xa (correspondant à y, qui, par suite de l'état 
variable de Tinstrument ou pour une cause analogue, seraient inévitablement les mêmes 
pour chaque observation faite avec cet instrument dans la temps ta^ mais qui varient 
accidentellement d'un moment à l'autre, de sorte qu'on a les équations 

aa — ^«+i'= (2) 
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avec les poids correspondants v, comme plus haut; les autres qui, accidentellement et 

indépendamment les unes des autres, influent sur les différentes observations de manière 
que leur action combinée donne pour les valeurs véritablement observées z^ au lieu de C«9 
et par suite les équations 

Za = ^a+pfa + qga (18) 

avec les poids Va correspondant aux équations (3) dans le problème sous sa forme plus 
simple. Avec n-|- 1 observations (indices de à n) on disposera donc de 2n-{- l équa- 
tions pour la détermination de n-f-S inconnues ^oi ^i.. .â?,, p et q. D*après la méthode 
des moindres carrés on aura alors 

ro^o = (vo + »)«o— va!i -f- + vofop + v^g^q 

0.1 0,1 

n^t = _ t7 «0+ («? + Vi + V)Xi—V a?2+ + VifiP + ViQiq 

0,1 0,1 1,2 1,2 



(19) 



Vo^a = —V œ^^i + (V + 1?.+ V) Xa — V ^«^.l+ . . . + V^aP + M«? 

a— 1,« a— l,o a,a+l «»«+! 

Vn^n = —V ^,-1 + (» + Vn) X^ + V»/„p + «^nfi^»? 

«— 1,« »*-l,n 

[t?/;t] = ro/o^o + t?i/l^i + + t?«/« A-« + + V^fnOin + (t?/*]i> + \^fg\ q 

\vgz\ = roJ^o^oH- Vl^'l^l + + «?aS'a«a+ + ^«.^««i» + [«/^IjP + I^*J? 

oii en général [t;r«] = t?oro«, + riri«i-|- • • . • + ^n*'»i^«* ^^t^^ Télimination entre ces 
équations nous employons toujours les valeurs de t/. et de u déterminées par (6) et (7), 

a,«-hl 

de sorte que la fraction continue (10) joue également ici un rôle essentiel. Nous nous 
servons en outre des mêmes quantités auxiliaires qui ont été désignées plus haut par y, 
mais que nous désignerons ici autrement, à savoir par \z\\ on aura ainsi \z^\ = zo et 

Va V 

«.«+1 

plus deux couples de valeurs {f] et [g] déduites d'une manière analogue de/ et de g et 
déterminées par (/oj =/o, [go] = go et par 

Ua V 

«t«+i 



Ua-f-V 
«,a+l 

«,«+1 



23 



et en outre quelques quantilés auxiliaires plus compliquées, quoique aos 

(/2(> W«t, if], {f9i o« t^*i 
toutes les cinq calculées d'après la même formule générale récurrente 

«-+if.+i»«+i( = t>.^-l . r,+, . «.+, + u.ir^,} — {u,— u) { r.\ 
avec la valeur initiale ir^it] = n««- 

Par l'emploi de ces quantités auxiliaires les n -^ I premières 
deviendront 






et doDDcroDt, pour la détermination de p et de q, les deux équations 

«.ll/,e,l - (Ah«-i1 = «.((/-*( — l/.t*)P + «.(l/^«! — i/-! 
M.((ff-M — lWt«»i) = "-(t/.?-! — (/-US-))'' + ""('?-*' ""U 

qui sont analogues aux équations qu'on aurait trouvées par la méthode 
fap-{-ga<i avait été une simple errear accidentelle, par conséquent à 

\vfz\ = [vP\p + \vfg\q 
\vgz\ = [vfg\p + [v9'^\q. 

On peut donner aux formules récuyrrentes pour les valenrs auxiliaires 
simples. Pour {r.^} ou posera, conformément à {8) 



et pour fd^-i^a-)-!) on a 

w.-h(|''"+i*-+i1 — V-\ l«-i} = frt-ifr.+iJ'.+i — ir.i \s,\) + u.[{r^,\ - 
ou 

«■+i(l»-.+.M-il — '■-+i«--n) + w ('■-+i*.+i — C.) I».t) + «-(''■«) '*■' " 



.••-t^- ■. 'r-* 
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et, en introduisant une nouvelle série de quantités auxiliaires, on peut en tirer des formules 
qui se prêtent à un calcul pratique, de sorte que (RS) est déterminé par 

a 
U ((RS) — {Va} [8a]) = Ua[{r^a) — (»•«) {8a)) (23) 

tandis qu'à Taide de cette équation, on calcule {ra+iSa^) par 

Wa+l(l^a+l»o+l) — ^a+l«a+l) = « ((K5) — ra+i««i-l) (24) 

a,a-i-l a 

Nous donnerons encore la formule 

qui peut avoir de l'importance comme contrôle pour (23)^ (24) et (22). 

Après avoir établi les formules fondamentales du problème, je ferai seulement 
remarquer que si fa est constamment =/ pour tous les indices, on aura aussi {/«J =/ 

et Ifa^a] =/{''«}) <i'oû il suit que le terme constant,^, est indissolublement lié aux a, 
en sorte que ^a+^ prend la place de au dans les équations (20) dont ravant-dernier 
terme disparaît, tandis que la première des équations (21) devient une identité et la 
seconde se simpliQe. Par conséquent, si Téquation du phénomène observé renferme un 
terme constant, celui-ci ne pourra être séparé des erreurs a. Cette solution s'applique 
donc sans changement aux équations de la forme Ca = o+p/a^ qga> 

Exemple de calcul. 

Il ne sera sans doute pas inutile d'éclaircir par un exemple le cas le plus simple 
de toute cette théorie. A cet effet je choisis une série d'observations ayant pour but de 
déterminer l'erreur de l'index du micromètre du cercle méridien de Copenhague, et qui 
ont été faites à Toccasion de la détermination, pendant l'été de 1879, de la différence 
de longitude entre Lund et Copenhague. J'ai attribué le même poids à chaque observa- 
tion de la coïncidence entre le fil mobile et le fil central fixe; je suppose en outre que le 
point de coïncidence entre les deux fils, lequel évidemment n'a pas été complètement 
invariable, a varié d'un moment à l'autre d'une manière purement accidentelle, et que sa 
variation moyenne a été la même dans chaque élément de temps, et par conséquent t; 

inversement proportionnel à ^a+i — ^at que le micromètre ait été employé on non dans 
l'intervalle; mais, comme on le verra, jai essayé de vérifier cette dernière hypothèse par 
une recherche spéciale. Je suis naturellement parti de la double supposition que les deux 
lois des erreurs étaient exponentielles et exemptes d'erreurs systématiques, mais il est 
douteux que la dernière ait été complètement justifiée. 

Comme unité pour les observations j'ai pris Vioo de tour de vis = 0*2633, et, 
pour économiser la place, j'ai rejeté les dixaines, centaines etc. de cette unité qui sont 
communes à toutes les observations. D'après cela, les unités de poids Ep et U^ doivent 
être déterminées de manière que Ep = 1 signifie que l'erreur moyenne d'une coïncidence 
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était 0*2633 , tandis que £« = 1 exige que le déplacement moyen de Tinstrument dans un 

jour ait été de 0*2633. Ep a provisoirement été déterminé par les répétitions directes 
réduites en moyennes, comme on a 36 équations indépendantes les unes des autres, cha- 
cune avec le même poids qu'une coïncidence, qui ont fourni une somme de carrés des 
écarts = 2*902, ce qui a donné E^ = 12*4. 

Par une compensation provisoire qui reposait sur des suppositions en partie in- 
exactes, j'avais trouvé Ep = 12*5 et -B,, = 24*8. Avec ces valeurs j'ai calculé celles de 

Va et de v, qui avec Za constituent les données dans la suite du calcul. Pour simplifier 
«,«+1 

ce calcul, j'ai en réalité opéré partout avec les doubles poids., mais pour plus de clarté 
je donne les résultats du calcul comme si cela n'avait pas eu lieu. Une exactitude de 2 
chifires a été suffisante. 
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Si la somme ^v^iPM + ^ v Af (a?,— a?«t*l °'®8t pas exactement = 74, il 

faut surtout Tattribuer à la détermination incertaine des quantités entre parenthèses dans 
le second terme de la somme, lequel correspond à des intervalles très petits. Pour la 
détermination de E'^ on a donc, dans 74 — 17*77 = 56*23 équations en eicès, une somme 
de carrés = 53'22, et en outre, dans 36 équations, une somme de carrés = I2'5 x 2*902 
= 36*27 qui provient de la réduction des observations en moyennes, soit en tout 92*23 
équations avec une somme de carrés = 89*46, ce qui donne E'^ = 1*03. Pour J?'», par 
contre, on a 73 — 56*23 = 16*77 équations en excès et une somme de carrés = 16*73, 
et par conséquent E'^ «= rOO. L*accord est donc satisfaisant, et je suis arrivé à ce 
résultat que le poids absolu d'une observation de coïncidence est 12*9, correspondant à 

une erreur moyenne = + 0*278 = + 0073, tandis que la stabilité du micromètre dans 
1 jour sidéral se mesure avec le poids absolu 23*8, correspondant à un écart moyen dans 

le même temps = +0*205 = +0*054. Relativement à notre supposition que les poids 
doivent simplement avoir été inversement proportionnels à la longueur des intervalles, on 
voit par le tableau suivant 

Intervalle. Equations en excès. Sommes des carrés. Ew. 

plus petit que 0*5 jour. . . 3*69 4*91 0*75 

entre 0*5 et 1 jour 7*08 615 115 

plus grand que 1 jour . . . 6*00 5*67 1*06 
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que les poids semblent bien devoir être un peu plus faibles dans les courts intervalles où 
le micromètre a été très employé, mais que ce résultat, par suite du petit nombre d'équa- 
tions en excès, est si incertain qu'il n'est nullement nécessaire de modiQer la supposition 
que Tusage plus ou moins grand qu'on a fait du micromètre a été sans influence sur sa 
stabilité. La question de savoir si, après cette compensation, il eiiste des traces d'erreurs 
systématiques ne peut non plus, pour le même motif, être définitivement répondue. Le 
phénomène qui, dans les tableaux des erreurs, semblerait pouvoir indiquer qu'elles 
s'écartent de la loi exponentielle ,* à savoir le petit nombre des variations de signes dans 
les différences Xa — ^^+i , ne saurait guère prouver un pareil écart ; mais, d'un autre côté, 
il est bien vraisemblable que les variations du micromètre peuvent en majeure partie être 
dues à celles de la température, et que le phénomène en question peut élre en connexion 
avec ces dernières. Toutefois, eût-on à sa disposition les observations thermométriques 
nécessaires, il ne vaudrait guère la peine d'y avoir égard. 

Pour éclaircir par un exemple l'usage de l'équation (16), je calculerai encore le 
point de coïncidence et son erreur moyenne à posteriori pour le long intervalle entre le 
21 et le 26 juillet. On trouve ici pour les temps intermédiaires U 

26-66 — ti 
P = 4-87 

«*== 2-51/> + 3-29(l— p) 
= 3-29 - O'ISp 
i/5(a?^) = 0-0208p + 0-0U2(l^/>) + 0-l724p(l — p) 
= 00142 + Ol790p - 01724p«, 

et par conséquent les déterminations suivantes: 

i X Limites probables 

21-79 2-51 ±0-14 2-37 à 2-65 

2200 2-54 ±017 237 - 271 

2300 2-70 ±0-23 2-47 - 293 

2400 2-86 ±0-25 2-61 - 3-11 

25-00 302+0-23 2-79 - 3-25 

26-00 3-18 ±0- 19 2-99 - 337 

26-66 3-29 + 0-12 3-17 - 3-41. 

L'incertitude atteint son maximum le 24*13 juillet avec une erreur moyenne de 
± 0*246 , et la détermination qu'on en tire pour x dans cet intervalle est donc partout 
plus exacte que celle qu'on aurait obtenue par une observation isolée de la coïncidence. 



